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Modifications of the infinite square well E(5) and X(5) descriptions of transitional nuclear structure 
are considered. The eigenproblem for a potential with linear sloped walls is solved. The consequences 
of the introduction of sloped walls and of a quadratic transition operator are investigated. 
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I. INTRODUCTION 

The E(5) and X(5) models have been proposed by 
Iachello 0> Q] to describe the essential characteristics of 
shape-transitional forms of quadrupole collective struc- 
ture in nuclei. The E(5) model, for 7-soft nuclei, and 
the X(5) model, for axially symmetric nuclei, are both 
based upon the approximation of the potential energy 
as a square well in the Bohr deformation variable /3. 
These models produce predictions for level energy spac- 
ings and electromagnetic transition strengths intermedi- 
ate between those for spherical oscillator structure and 
for deformed 7-soft [j| or deformed axially-symmetric ro- 
tor structures. 

The X(5) predictions for level energy spacings and elec- 
tromagnetic transition strengths have been extensively 
compared with data for nuclei in transitional regions be- 
tween spherical and rotor structure [U IS E IS IS fTot 
ITll IT2I. fl^ l . For several such nuclei, including the 7V=90 
isotopes of Nd, Sm, Gd, and Dy, the X(5) predictions 
match well the yrast band level energies and the excita- 
tion energy of the K*=0$ band head [Fig. [IJ^b)]. The 
X(5) predictions also reproduce essential features of the 
electric quadrupole transitions from the ^^=0^" band 
to the ground state band: the presence of strong spin- 
ascending interband transitions but highly-suppressed 
spin-descending transitions. 

However, several discrepancies exist between the X(5) 
predictions and observed values. The spacing of level 
energies in the K^—O^ band is predicted to be much 
larger than in the ground state band, but empirically 
at most a slightly larger energy scale is found for the 
K*=Q% band [Fig. [Jc)] ISEllll This overprediction 
is encountered in descriptions of transitional nuclei with 
the interacting boson model (IBM) and geometric col- 
lective model (GCM) as well [H For nuclei with 
yrast band level energies matching the X(5) predictions, 
the yrast band B(E2) strengths tend to fall below the 
X(5) predictions, and sometimes even below the pure ro- 
tor predictions (see Fig. 2 of Ref. [TJ). For the A=90 
nuclei, the transitions between the if 7r =0^" and ground 
state bands have strength ratios typically matching those 
predicted, but their strength scale is considerably weaker 
than predicted H @, Hfl ^ . 

It is thus necessary to ascertain which aspects of the 



X(5) description are most important in determining the 
predictions for these basic observables. The square well 
potential involves an infinitcly-stcep "wall" in the poten- 
tial as a function of [3, presumably a radical approxima- 
tion. Moreover, the model has so far been used only 
with a first-order electric quadrupole transition oper- 
ator, but the likely importance of second-order effects 
has been noted by Arias |22j and by Pietralla and Gor- 
bachenko J2j|. In the present work, the infinitely stiff 
confining wall is replaced with a gentler, sloped wall, con- 
structed using a linear potential. The effects upon cal- 
culated observables of the introduction of a sloped wall 
and of a quadratic transition operator are addressed. A 
computer code for solution of the sloped well eigenprob- 
lem is provided through the Electronic Physics Auxiliary 
Publication Service |24|. 



II. SOLUTION METHOD 

Consider the Bohr Hamiltonian £| 
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where [3 and 7 are the Bohr deformation variables and 
the M K are angular momentum operators, with potential 
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Since this potential is a function of f3 only, the five- 
dimensional analogue of the central force problem arises. 
The usual separation of "radial" (f3) and "angular" vari- 
ables |S [23 occurs, yielding eigenfunctions of the form 
*(/3,7,w)=/(/3)$(7,w), where w=(0i, i? 2 , #3) are the 
Euler angles. The angular wave functions ^(7, uj), com- 
mon to all 7-independent problems, are known [26j. For 
the radial problem, following Rakavy 25], it is most con- 
venient to work with the "auxiliary" radial wave function 
(p(/3)=(3 2 f((3). This function obeys a one-dimensional 
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FIG. 1: Evolution of the (a) 4+ energy, (b) Oj energy, and 
(c) energy spacing scale of the excited + sequence, all nor- 
malized to the 2J" energy, across the iV=90 transition region, 
for the Nd (□), Sm (O), Gd (A), and Dy (♦) isotopic chains. 
Shown for comparison are the X(5) predictions (solid line) 
and the present sloped well predictions for various values of 
the parameter S defined in 11211 — S— 100 (long-dashed line), 
S=50 (short-dashed line), and 5*=25 (dotted line). Data are 
from Refs. El [H [H Ea . (Figure based upon Ref. 0.) 



Schrodinger equation with a "centrifugal" term 
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where the centrifugal coefficient a is related to the 
0(5) separation constant r (r=0, 1,...) by a=(r + 



1)(t + 2). For problems with a more general potential 
V(J3,i)=Vp(fi) + Vjfr), Iachello showed that an ap- 
proximate separation of variables occurs, provided that 
V-fi^y) confines the nucleus to 7«0 (see Ref. 2j for de- 
tails). In this "7-stabilized" case, the eigenfunctions 
are of the form ^((3, 7, w)oc/(/3)?7(7)0klm(w), where the 
4>klm{w) are the conventional rigid rotor angular wave 
functions Q for angular momentum L, z-axis projec- 
tion M, and symmetry axis projection K . The auxil- 
iary radial wave function again obeys 0, but now with 
a=\L{L + 1) + 2. 

In the region (3<{3 W , the potential V((3) of @ van- 
ishes, and the radial equation reduces to the Bessel 
equation of order v=(a + 1/4) 1 / 2 . The solutions with 
the correct convergence properties at the origin are 
^(/3)oc^ 1 /2j [/ ( £ i/2 ; g) ) w here e={2B/h 2 )E. In the region 
/3>/3 w , where the potential is linear in /3, an analytic so- 
lution does not exist for the full problem with centrifugal 
term. For a=0 only, reduces to the Airy equation, 
with solutions y(/3)ocAi[c 1 / 3 (/3 — (3 W ) — c~ 2 / 3 e], where 
c={2B/h 2 )C. 

The analytic solutions obtained for a=0 provide a very 
efficient basis for numerical diagonalization to obtain the 
true a^O solutions of the radial equation J2J. It is first 
necessary to obtain a basis set of a—0 solutions 
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The eigenvalues of e are determined by the condition that 
(p(/3) be continuous and smooth at the matching point 
f3=f3 w . This yields a transcendental equation which is 
solved numerically for s. The normalization coefficients 
Ni and N2 then follow from continuity and the require- 
ment J °° dj3\ip(l3)\ 2 =l. Since the radial equation J5J has 
the form of a one-dimensional Schrodinger equation, its 
solution for general values of a may be carried out as 
the matrix diagonalization problem for a corresponding 
"Hamiltonian" matrix h, including the centrifugal po- 
tential, with respect to these a=0 basis functions, with 
entries 
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d/3<pr (P)^<pr°w- (5) 

Convergence in this basis is rapid — for instance, the 
eigenvalues of the ground state and first excited radial 
solution converge to within ~1.5% of their true values 
with a truncated basis of only 5 eigenfunctions. Values 
shown in this paper are calculated for a basis size of 25. 
For illustration, an example potential, with centrifugal 
contribution, and the corresponding calculated eigenval- 
ues are shown in Fig. [21 

Electromagnetic transition strengths can be calculated 
from the matrix elements of the collective multipole op- 
erators. The general E2 operator for the geometric 
model [27l I2H. I23 may be expanded in laboratory frame 
as M 
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FIG. 2: Energies of low-lying + and 2 + levels for the sloped 
well potential with S—50. The potential without the five- 
dimensional centrifugal term is shown (solid curve), together 
with the potential including the centrifugal contributions for 
L=0 and L—2 (dashed curves). The energies of the corre- 
sponding states for the X(5) model are shown for comparison 
at right. [The S=50 calculation is for 2B/fi 2 — l and j3 w =l, 
while the X(5) calculation is scaled to p m — lA0 to provide the 
same ground state eigenvalue.] 



For the present purposes, it is necessary to reexpress this 
operator in terms of the intrinsic frame coordinates and 
D 2 (co) [3l|, giving, to second order in /3, 

M(E2; ,i) = Ai/3[DlS cos 7 + fypfc + K-*) ^7] 

- ^1 A 2 f [Dl* cos 2 7 - \{Dl$ + Dl*_ 2 ) sin 2 7 ] . 

(7) 

In both the 7-independent and 7-stabilized cases, the 
matrix element of VJl(E2; fi) between two eigenstates fac- 
tors into an angular integral and a radial integral. Here 
we consider matrix elements between unsymmetrized 7- 
stabilized wave functions 41 
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Dl; K (cu)<i> aK j((3 n ), (8) 



as needed in calculations for the rigid rotor, X(5), or 7- 
stabilized sloped well models. The matrix element sepa- 
rates into intrinsic and Euler angle integrals, yielding 



Aih - \I-A2I2 



{^a'K'J'M' I \M(E2; fi)\ \^ a RJM 
(-) J '- J (J'K'2(K - K')\JK) 

in terms of 

h = J dr$;, A -, J ,(/3,7)/3cos 7 $ Q Kj(/3,7) 

h = Jdr K>K>AP>l) P 2 cos27$ QXJ (/3, 7 ) 



(9) 
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for K'-K=0 or 



/, = / dr$^, J ,(/3,7)-^=/3sin 7 $ Q/fJ (/3,7) 

dr K'k>j> (A 7) (-)^ 2 sin 2 7 *aKj(P, 7), 

(11) 

for K'— K=±2, where dr=P 4 d(3 |sin 3 7 | d 7 and the 
reduced matrix clement normalization convention 
is that of Rose [32 • (The matrix elements of the 
symmetrized wave functions, for Kj^O, may be 
calculated from this matrix element as usual 

Bo 

Considering the present /3- 7 separated wave functions 
$(/3, 7 )=/(/3)?7( 7 ), for the case of no 7 excitation (so 
K'=K=0), and under the approximation 7~0, these 
integrals reduce to h=j l3 i dj3f a > K 'j'{l3)!3f a Kj{P) 
and h=J f3 4 df3f a , K ,,j,({3)[3 2 f aK j(0). Transition 
strengths are B(E2; J^J')=(2J' + 1)/(2J + 1) 
I (J'||SDt(.E2)||J) | 2 . Quadrupole moments, defined 
by e g J =(167r/5) 1 / 2 (JJ\Tl(E2; 0)| J J), may be calcu- 
lated as eQ J =(167r/5) 1 / 2 (JJ20|JJ) (J\\Wl(E2)\\ J). 

The following calculations can be considerably simpli- 
fied if it is noted that the eigenvalue spectrum and wave 
functions depend upon the Hamiltonian parameters B, 
(3 W , and C only in the combination 



S — ^2 ftwC-i 



(12) 



to within an overall normalization factor on the eigenval- 
ues and overall dilation of all wave functions with respect 
to (3. [This follows from invariance of the Schrodinger 
equation solutions under multiplication of the Hamilto- 
nian by a constant factor and under a transformation of 
the potential V {(3)=a 2 V{a(i) [33|.] For a given value of 
S, the numerical solution need only be obtained once, 
at some "reference" choice of parameters (e.g., 2B/h 2 =\ 
and p w — l), and the solution for any other well of the 
same S can be deduced analytically. Specifically, sup- 
pose the reference calculation yields an eigenvalue e and 
a normalized radial wave function /(/?). Then a calcula- 
tion performed for the same B and S but for a different 
width P' w produces the eigenvalue e' and normalized wave 
function /'(/3) given by the simple rescalings 



e> = 0- 2 e 
/'(/?) = f3'~ 5 / 2 f(P/f3' w ), 



(13) 



and the radial integrals scale to /{ = (3' w Ii and I' 2 = 
Thus, the essential parameter which controls the relative 
strengths of the linear and quadratic terms of the E2 
operator is A*=Ai{i w IA\, in terms of which the matrix 
element in (0 is 



) J '- J (J'K'2(K - K')\JK) 
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(14) 
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Ratios of E2 matrix elements depend only upon S and 
A'. 

A computer code for solution of the sloped well eigen- 
problem and for calculation of the radial matrix elements 
between eigenstates is provided through the Electronic 
Physics Auxiliary Publication Service j^. This code 
also calculates observables for the E(5) and X(5) models. 



III. RESULTS 

In the following discussion, let us restrict our atten- 
tion to 7-stabilized structure relatively close to the X(5) 
limit of the sloped well model, since this regime is most 
directly relevant to the transitional nuclei recently con- 
sidered in the context of the X(5) model. The sloped 
well potential approaches a pure linear potential as [3 W 
vanishes at fixed slope (that is, as S— >0) and approaches 
a square well as the slope goes to infinity at fixed j3 w 
(that is, as S— >oo). It can thus produce a much wider 
variety of structures than are considered in the present 
discussion. However, calculations for the full range of 
these cases may be obtained with the provided computer 
code 0. 

First we examine the energy spectrum, comparing it 
to the X(5) spectrum. Naturally, the eigenvalues for the 
sloped well are lowered relative to those for the X(5) well 
of the same (3 W , as the outward slope of the wall effec- 
tively widens the well, causing level energies to "settle" 
lower. The essential feature is that the widening of the 
well introduced by the wall slope is a relatively small 
fraction of the well width at low energies, while it is 
much greater at high energies, as may be seen by in- 
spection of the potential (Fig. |2J). Thus, the high- lying 
levels experience a disproportionately greater increase in 
the accessible range of /3-values than do low-lying levels 
and consequently are lowered in energy relative to the 
low-lying levels. 

From the calculated energies, it is seen that as S is de- 
creased from infinity the higher-spin levels within a band 
are lowered more rapidly than the lower-spin members, 
resulting in a reduction of the ratio Rn/2=E{^[) / E{2^[ ) 
for the yrast band [Fig.QJa)] and a lowering of the curve 
of E versus J for each band (Fig. |3J) . The excited band 
head energies are lowered as well [Fig. db)]. But the 
most dramatic change is the rapid collapse of the spacing 
scale of levels within the excited bands relative to that of 
the ground state band [Figs.^c) and|3|. For 5«50, the 
predicted energy spacing scale within the _R" n '=0j band 
is reduced sufficiently to be consistent with the spacings 
found for the 7V=90 transitional nuclei, while the ener- 
gies of low-spin yrast band members and the ^=0^ 
band head are still relatively close to their X(5) values, 
as shown in Fig. |3J 

The second order term in the E2 operator J7J) can in- 
terfere either constructively or destructively with the first 
order term. For all transitions between low-lying levels 
considered here, the radial integrals I\ and I2 in l|14|) 
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FIG. 3: Yrast and K n = band level energies, normalized 
to E(2f), for 150 Nd (□), 152 Sm (O), 154 Gd (A), 156 Dy (♦), 
and 162 Yb (O). The predictions for X(5) (solid curve) and 
the sloped well with 5=50 (dashed curve) are shown for com- 
parison. Data are from Refs. El EI El El El 



have the same sign. Thus, negative values of A' lead 
to constructive interference [note the negative coefficient 
in (|14fl ] , while positive values lead to destructive interfer- 
ence. For the X(5) square well, the higher-spin members 
of the yrast band have larger average (3 values than do 
the low-spin members, so the quadratic term is relatively 
more important for the higher-spin levels. In the case of 
destructive interference, the curve showing the spin de- 
pendence of B{E2) values, normalized to B(E2; 2\ — > 
0+ ), falls below that obtained with the simple linear E2 
operator, as seen in Fig. Ufa). The broad r ang e of such 
curves obtained experimentally (see Ref. [TlJ ) can be 
qualitatively reproduced with different values of A'. De- 
structive interference also reduces the interband B(E2) 
strengths and the in-band B{E2) strengths within the 
K*=0% band, relative to B(E2;2+ -> 0+) [Fig. EJb)] , 
ameliorating the overprediction of interband strengths in 
the X(5) model. The spin-descending interband tran- 
sitions in the X(5) model have highly-suppressed linear 
E2 matrix elements, so these transitions are very sensi- 
tive to even a small quadratic contribution. Values of 
A' which give only moderate modifications to the other 
transitions can give complete destructive interference for 
these spin-descending transitions. The spin dependence 
of quadrupolc moments within the yrast band is shown 
in Fig.|S». 

Observe that the situation just described differs consid- 
erably from that encountered for a pure rotor. For a rigid 
rotor, the intrinsic wave function ^ a x(P, 7) is the same 
for all levels within a band, so 1% provides only a uni- 
form adjustment to the intrinsic matrix element between 
bands. Inclusion of the second order term in 9Jt(_E2) 
thus leaves unchanged the ratio of any two B(E2) val- 
ues within a band or the ratio of any two B(E2) values 
between the same two bands. 

Although inclusion of a quadratic term in the E2 oper- 
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FIG. 4: Yrast band B(E2) strengths and squared quadrupole moments, normalized to B(E2;2f — ► Of). (a,b) Values for the 
X(5) model (S^oo) calculated with a quadratic E2 transition operator, for A' ranging from —1 to 2 in equal steps. (c,d) Values 
for the sloped well (.5=50, 100, 200, 500, and oo) calculated with a linear E2 transition operator. 
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FIG. 5: Level schemes and selected B(E2) strengths for (a) the X(5) model with a linear E2 operator, (b) the X(5) model 
with a quadratic E2 operator (A'=0.7), and (c) the sloped well (S'=50) with a linear E2 operator. Arrow thicknesses are 
proportional to the logarithm of the B(E2) strength. 
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FIG. 6: Level scheme and selected B{E2) strengths (a) for 
the sloped well with parameters chosen to approximately re- 
produce the observed low-energy structure of 150 Nd (5=75, 
A' =0.6) and (b) as measured for 150 Nd 0, G3- Arrow 
thicknesses are proportional to the logarithm of the B(E2) 
strength. Limits are indicated on experimental B(E2) 
strengths for transitions with unknown E2/M1 mixing ratios. 



ator with A'>0 can at least qualitatively explain the dis- 
crepancies between the X(5) B(E2) predictions and em- 
pirical values, this explanation is not entirely satisfactory. 
Many different spin dependences of the B{E2) values 
within the yrast band are observed for nuclei with simi- 
lar energy spectra , and these require correspondingly 
varied, apparently ad hoc choices of the parameter A 1 for 
their reproduction. Moreover, it is possible to obtain esti- 
mates for the coefficients A\ and Ai in the geometric E2 
operator based on a simple model of the nuclear charge 
and current distribution, as described in Refs. [28l |29|. 
and these values yield A' fa— 0.2, giving weak constructive 
interference for the low- lying transitions. In the interact- 
ing boson model, the E2 transition operator is of the 
form T( B2 >oc(d t x s + s 1 * x J) (2) +x(d t x d)( 2 \ in terms of 
the boson creation operators and d) , where the value 
X = ~ v7 /2 is commonly used in calculations involving the 
transition from spherical to axially-symmetric deformed 
structure [35|,|3(|. In the classical limit, {S xs+stx d)^ 
may be approximately identified with the linear term of 
the geometric model transition operator and x(d^ x d)^ 
with the quadratic term. The addition of these terms 
is constructive for low-lying transitions, and the relative 
contribution of the second-order term is comparable to 
that obtained for A'pz— 1 in the present description. 

The effect of sloped walls on the calculated B(E2) 
strengths is dominated by the greater broadening of the 
well at high energies than at low energies discussed above. 
While all the eigenfunctions "spread" in f3 extent rela- 
tive to those for the square well, this spreading is most 
pronounced for the high- lying levels. Since the first or- 
der E2 operator is proportional to /3, the E2 matrix ele- 
ments tend to be enhanced for the higher-lying levels. In 
the yrast band, the in-band B(E2) strengths for higher- 



spin band members are increased relative to those for 
the lower-spin band members, as are the quadrupole mo- 
ments for higher-spin band members [Fig. EIc,d)]. Sev- 
eral of the interband B(E2) strengths are also increased 
relative to B(E2;2^ -> Of) (Fig. 0). The changes in 
B(E2) values induced by decreasing 5* are largely oppo- 
site in sense to those produced by introduction of the 
second-order term in the E2 operator. The parameters 
S and A' may be chosen so as to balance these two effects 
against each other, except that for the spin-descending 
interband transitions the strong destructive interference 
tends to dominate. 

To allow comparison with empirical values, in Fig.^a) 
predictions obtained with the sloped wall potential and 
quadratic E2 operator are shown for parameter values 
chosen to approximately reproduce the observed low- 
energy structure of 150 Nd. The experimental values are 
given in Fig.HJb). 

Finally, let us consider the effects of wall slope on 
the properties of the K 7r =2^ band, or 7 band. Within 
the 7-stabilized separation of variables of Ref. 0, the 
properties of this band are largely independent of the 
specific choice of 7-confining potential ^(7). This po- 
tential determines the band head energy as well as 
the 7-dependent wave function 77(7). The wave func- 
tion, however, simply contributes a normalization factor 
J I sin 37 1 £^7771 (7) sin 7770(7) to l\ in (|llfl . and an analo- 
gous factor to I2, common to all electromagnetic matrix 
elements between the K v =2f band and the K n ~0^ and 
O2 bands. Although these quantities can be calculated 
for any particular hypothesized form for Vy(j), such as 
a harmonic oscillator potential Q, they in practice may 
be treated as free parameters. 

The essential feature of the K 1T =2f band is that the ra- 
dial wave function for each of its members is the "ground 
state" solution of the radial equation J3J for the given 
angular momentum. This K^=2\ band is thus essen- 
tially a duplicate of the yrast band, displaced to a higher 
energy by the excitation energy in the 7 degree of free- 
dom, with energy spacings and radial wave functions for 
the even spin members identical to those for the yrast 
band, but with the addition of odd spin members and 
with different angular wave functions. (Note that for 
K 7^0 Bijker et al. use a different separation proce- 
dure from that in Refs. @,[6j, yielding a modified form of 
the radial equation with a=^[L(L + 1) — K 2 ] + 2, which 
changes the energy spacings and in-band radial matrix 
elements by <5% relative to those of the yrast band.) 
Thus, the dependence of K 7r =2f band properties upon 
wall slope closely matches that of the yrast band prop- 
erties. Notably, the K 7r =2f band does not demonstrate 
the rapid decrease in energy spacing scale with decreasing 
wall slope exhibited by the ^^=0^" band, as illustrated 
in Fig.[7Ja). This is at least qualitatively consistent with 
the observed similarity of the yrast and K 71 ^^ , but not 
K n =(>2 , band energy spacings in the N=90 X(5) candi- 
date nuclei. 

Since the even spin members of the K 7V =2'^ band pos- 
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J J 

FIG. 7: Dependence of properties of the K V =2 J [ band upon wall slope. Values are shown for the X(5) model (solid line), 
S^lOO (long-dashed line), 5=50 (short-dashed line), and S—25 (dotted line), (a) Yrast, K 71 =0^ , and K 1 * =2\ band energies, 
illustrating the identical dependences of the yrast and K n =2f bands. The K 71 T =2J band head energy is arbitrary (see text), 
(b) B(E2) branching ratios for AJ=±2 transitions from the K w =2^ band to the yrast band. (c,d) B(E2) strengths for 
AJ=±2 transitions from the K n —2~l band to the K 7 '=Q~\ band, normalized to B(E2;2^ — » 0^). The separation of variables 
of Refs. 0, has been used for all calculations. 



sess the same radial wave functions as the yrast band 
members, the strengths of transitions within the K 7T =2\ 
band or between this band and the yrast band depend 
upon the same radial matrix elements as do the yrast in- 
band transition strengths and quadrupole moments al- 
ready considered. Consequently, decreasing wall slope 
leads to a moderate enhancement of the interband tran- 
sition strengths involving higher spin levels, directly com- 
mensurate with the increases shown in Fig.0Jc,d). The 
dependence of branching ratios from the K 7T =2~[ band to 
the yrast band on wall slope is shown in Fig. EJb) . Tran- 
sitions between the K 71 =2^ band and the if^Oj band 
depend instead upon radial matrix elements which con- 
tribute to the X 7r =0^ to yrast band transition strengths. 
(The small radial matrix element values which yield the 
characteristic suppression of spin- descending transitions 
from the K n =02 band to the yrast band here yield 
a suppression of spin- ascending K^=2\ to 0^ transi- 
tions.) Decreasing wall slope yields enhancement of the 



allowed transitions and, for the low-spin levels, either 
little change or substantial reduction of the suppressed 
transitions [Fig. E{c,d)]. Detailed quantitative predic- 
tions for the K 7r =2~[ band level energies and electromag- 
netic observables, using either the separation of variables 
of Refs. or that of Ref. 0], ma y be obtained with 
the provided code |24j . 



IV. CONCLUSION 

The use of a /3-soft potential within the geometric pic- 
ture has recently received attention as providing a simple 
description of nuclei intermediate between spherical and 
rigidly deformed structure. From the present results, it is 
seen that the energy spacing scale of states within excited 
bands is highly sensitive to the stiffness of the well bound- 
ary wall. A potential for which the well width increases 
with energy can produce a more compact spacing scale 
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for excited states than is obtained with a pure square 
well, providing much closer agreement with the observed 
energy spectra for nuclei in the 7V«90 transition region. 
It is also found that a second-order contribution to the 
E2 transition operator can lead to a wide range of pos- 
sible yrast band B(E2) spin dependences, as well as to 
modifications of off-yrast matrix elements. However, a 
systematic understanding of the proper strength for this 
second-order contribution is needed if the E2 operator is 
to be applied effectively. 
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